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A selection of recent heavy-flavor results from OPAL using the LEP1 data sample are pre- 
sented. The average polarization of b baryons in hadronic Z° decay has been measured to be 
— 0.56to'i3 (stat.) ± 0.09(syst.) using semileptonic decays of Ab baryons. A search has been con- 
ducted for the radially excited D*' and has produced a 95% CL upper limit on its production 
of /(Z° -> D* /:t (2629)) x Br(D*' ± -> D^tt+O < 2 - 1 x w ~ 3 ■ Finally, the measurement of 
the product branching ratio ,f(b -> Ab) x Br(A b -> AX) = (2.67 ± 0.38(stat)+°;gp(syst.))% has 
been made. This measurement, along with an earlier measurement of the product branching ratio 
/(b — ■» Ab) x Br(Ab — > AiX), has been used to compute an updated Rm = Br(Ab — > A£X)/Br(Ab — » 
AX) = (8.0± 1.2(stat.) ±0.9(syst.))%, consistent with the expected low semileptonic branching frac- 
tion of the Ab inferred from its short lifetime compared to the other b hadrons. 



0\ 



I. MEASUREMENT OF THE AVERAGE POLARIZATION OF B BARYONS IN Z° DECAYS 



In the Standard Model, b quarks produced in e + e collisions with y/s « m z o have an average longitudinal polar- 
ization of —0.94; the sign indicates that the b quarks are predominantly left-handed. The variation of this polar- 
ization with the angle of emission of the b quark 9\> [] is only ±2%, and the dependence on sin 2 #yK is fairly weak: 
" <9(-Pb)/<9sin 2 6w ~ 0.63. This average polarization of b quarks is the quantity Ay, measured by SLD Q), and is an 
(f) ■ ingredient to the A^ B measurements from LEP || . 

In the heavy-quark limit, the spin of the b quark is expected to decouple from the light degrees of freedom produced 
by the fragmentation processes. The process of stable hadron formation has a profound effect on the polarization 
0^ 1 retained in the final state. The pseudoscalar mesons B° and carry no polarization information. The vector meson 
B* decays by photon radiation which couples mainly to the light degrees of freedom within the meson, and hence 
the radiated photons are very nearly isotropic When directly produced, the Ab is expected to retain nearly the 
full polarization of the parent b quark Q , where the light valence quarks form a spin-zero system. Some Ab mesons 
are produced from strong decays of the Sb and the ££. The energy difference between these states is the hyperfine 
*£h 1 splitting induced by the relative orientation of the b quark spin and the spin-1 diquark. These states are produced in 
, a coherent superposition in the hadronization process. If the Sb and the have narrower widths than their mass 
splittings, then they decay incoherently; their lifetime is long enough to flip the spin of the b quark. If the states 
5_i ' are broader than the splitting, then the b polarization retention is higher. The widths of these states are predicted 
to be of the order of the mass splitting between them, hence substantial depolarization is expected for these cascade 
decays A measurement of the average polarization of Ab baryons tests both Heavy Quark Effective Theory 

(HQET) and also models of baryon formation. 

The polarization of Ab baryons is most easily measured using their semileptonic decays Ab — > i^V£K. To 
conserve angular momentum in the maximally parity- violating decay, the charged leptons l~ are preferentially emitted 
antiparallel to the spin of the Ab, while the TVs are preferentially emitted parallel to the spin. When boosted into the 
laboratory frame, the charged lepton energy is on average higher than the neutrino energy. The spectra depend on the 
Ab momentum spectrum and the details of the semileptonic decay, which are influenced by mb/m c , QCD corrections, 
and form factors. The ratio of the average charged lepton energy to the average neutrino energy partially cancels these 



x In the right-handed OPAL coordinate system, the positive z axis points along the e~ beam and the x axis points towards 
the center for the LEP ring. The polar and azimuthal angles are denoted by 6 and <j>, and the origin is taken to be the center 
of the detector. 

2 Charge conjugate processes are implied in this analysis and throughout this article. 
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FIG. 1. Invariant mass distribution of p7r~ combinations in the right-sign A£~ir + sample (open histogram) shown with the 
distribution for wrong-sign A£ + n + combinations (shaded histogram); the histograms are not added. The signal region is shown 
in brackets. The peak in the wrong-sign distribution arises from genuine lambdas from fragmentation or B meson decay paired 
with real or fake leptons. 




effects and is statistically more powerful than either energy alone |6j]. This analysis is based on the event-by-event 
distribution of the ratio E v /Ei. The ingredients needed for the analysis are a purified sample of semileptonic Ab 
decays, an estimation of the background contribution, and measurements of the lepton and neutrino energies. 

The OPAL measurement is described in detail in Reference and is summarized here. The selection uses correla- 
tions between reconstructed A — > p7r _ decays and high-(p,p t ) leptons similar to the selection of but with improved 
electron identification, different lepton kinematic cuts, and the inclusion of the requirement of a 7r + , whose charge is 
required to be correlated with the reconstructed A to take advantage of the decay chain Ab — > A+ — > An + X. The ir + 
is additionally required to be consistent with coming from the common A-lepton vertex. 

The invariant mass distribution of the selected p7r~ combinations is shown in Figure ^, separately for the A£~tt + 
(right-sign) and A£ + ir + (wrong-sign) combinations. From a total of 4.3 million hadronic Z° decays used in all three 
analyses presented in this article, a total of 912 right-sign A£ir + combinations are selected with an overall b baryon 
purity of 69%. A total of 316 wrong-sign combinations are also selected. The wrong-sign combinations were used to 
estimate the level of background as well as the shape of the background distribution of E v jEi. Part of the wrong-sign 
background consists of real leptons from Ab decay combined with fragmentation A's, which are preferentially produced 
with the wrong sign because of conservation of baryon number; these events contain the same polarization information 
as right-sign events. 

To estimate the neutrino energy, the missing energy in the signal hemisphere is computed as the difference between 
the expected total energy in the signal hemisphere and the measured energy in the same hemisphere: 

2 2 
III III 

rp Tp i hemi recoil rphemi /i \ 

E v - Ebeam H 77; E vis , (1) 

where Ebeam is the beam energy, E^f™ is the measured visible energy in the signal hemisphere, rrihemi is the invariant 
mass of the observed particles in the signal hemisphere and m reco u is the invariant mass of the opposite hemisphere. 
The estimated resolution of the neutrino energy is 3.5 GeV. The measurement error on the charged lepton momentum 
is negligible by comparison, although the momentum scale and spectrum modeling introduce a systematic uncertainty. 
The distribution of the measured E v is compared between data and Monte Carlo for events with a high-(p,p_L) lepton 
and also for those without, providing samples enriched in high-energy neutrinos and depleted in them. Differences in 
the E v spectra are used to estimate corrections and systematic errors. 
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FIG. 2. Bottom: distribution of the ratio of reconstructed energies E v jEg for the right-sign combinations (solid circles) and 
wrong-sign combinations (open circles). The solid and dashed curves are the best-fit models to the distributions in the right - 
and wrong-sign distributions, respectively. Top: fractional residuals in the fit to the right-sign sample, for the central value of 
—56%, and also showing models with (P^ h ) values of —1.0 and 0.0. 



TABLE I. Summary of systematic uncertainties in the measurement of {P^ h }. 



Source of Uncertainty 


A{P^) 


Source of Uncertainty 




E v resolution 


±0.02 


Ac polarization 


±0.02 


E v reconstruction 


±0.05 


b^ r 


±0.01 


Eg scale and shape 


±0.03 


Fitting method 


±0.03 


Selection criteria 


±0.02 


Theoretical uncertainty (form factor 


±0.03 


Background fraction and shape 


±0.04 


modeling, QCD corrections, m c /mb 




b fragmentation 


±0.03 







Total: ±0.09 



The binned fit to the measured E^/Ei distribution was performed by comparing it to fully simulated Monte Carlo 
samples generated with various values of the average polarization (P^ b ). The right-sign and wrong-sign spectra were 
fit simultaneously. These distributions and the fit are shown in Figure |^. The x 2 per degree of freedom is 0.91 and 
the fit value is (P^ b ) = — 0.561^13, where the errors are statistical. 

Table Q summarizes the sources of systematic uncertainty considered. Added together in quadrature, the total 
systematic error is estimated to be ±0.09. The average longitudinal polarization of b baryons is measured to be 

(p£») = -0.56+°;^(stat.) ± 0.09(syst.). (2) 

Including the effects of both statistical and systematic uncertainties, the 95% confidence level (CL) acceptance region 
is —0.13 > (Pi h ) > —0.87, excluding for the first time zero polarization at the 95% CL, and also disfavoring full 
preservation of the original polarization of —0.94 through the fragmentation and hadronization processes. 
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II. SEARCH FOR THE RADIALLY EXCITED 0*^(2629) 

In February 1998, the DELPHI collaboration reported evidence || for a radially excited charm meson, labeled the 
D*'. The lowest-energy radially excited state is called the D', but the D* /=t is easier to select because it is expected 
to decay to D* ± 7r + 7r _ with a significant branching ratio. DELPHI'S measured signal rate is expressed in terms of the 
rates of production of Dj 1 and Dj mesons in Z° decay: 

i? DEL PHi = (iVD " ±)Br(D * /+ ^ D * + ^ } = 0.48 ± 0.18(stat.) ± 0.10(syst.). (3) 
E J= x, 2 (^ D wo)Br(DW°^D*+7r-) 

The 95% CL upper limit on the width of the observed resonance is 15 MeV and is limited by the detector resolution. 
The observed mass is 2637±6 MeV. 

The expected mass of the D*' is 2629±20 MeV [|l0|, consistent with the DELPHI observation. An estimate of the 
decay width T(D*' — > D* + 7r + 7r~) using a harmonic oscillator approximation jll| is < 1 MeV, although a more recent 
computation by Melikhov and Pene |]l2|| of the same partial width using heavy-quark symmetry is between 120 and 
160 MeV. Melikhov and Pene also estimate the width T(D*' — > D^ir to be between 45 and 450 MeV, indicating some 
trouble in interpreting the observation. In addition, CLEO has searched for the D*' + in an analysis sensitive only to 
e+e" cc -^D*'+ |f] and has placed a limit of R < 0.16 at 90% CL. 

OPAL has searched for the D*' + and the selection procedure is summarized here. Charged D* + candidates are 
selected in the decay mode D* + -^D°7r + ,D° -^K~7r + . The K~7r + combination is required to have its invariant mass 
within ±2cr of the expected D° mass; the resolution is 25 MeV in both data and Monte Carlo. The D* + is formed 
with a pion of opposite sign to the kaon and the scaled energy of the D* + , xe — £b*+ /Ebeam is required to exceed 0.2. 
The mass difference Am = mo*+ — m^o must be between 142 and 149 MeV, corresponding to ±3<r in the resolution. 
Additional background rejection is obtained by placing requirements on the hclicity angle and the angle between the 
kaon and the D° flight direction in the D° rest frame. OPAL's excellent dE/dx capability |l4| is additionally used to 
purify the kaon sample. The D*' + candidates are formed by combining the D* + candidates with pairs of oppositely 
charged tracks, in order to explore DELPHI'S observation in the D* ± n + n^ decay mode. These pions have loose 
dE/dx requirements placed on them. The reconstructed mass of the D*' + is corrected by using the PDG value for 
the D* + mass and the measured difference between the D*' + and the D* + candidates. 

The analysis proceeds separately for events tagged as bottom and charm, in order to suppress background from 
Z° — » uu, dd, and ss events. To select candidates from a sample enriched in Z° — > cc events, the D* + candidate is 
required to have xe > 0.4, and the magnitude of the vector sum of the momenta of the two pions added to form the 
D*' + candidate is required to be greater than 3.6 GeV. To select candidates in Z° — > bb events, the separation of the 
D° vertex from the primary divided by its error is required to be larger than 4.0. A loose requirement is made on the 
separation between the D° vertex and the D*' + pion vertex to reduce the combinatoric contamination. Using a Monte 
Carlo simulation, the efficiency for selecting D*' + mesons in Z° — * cc events is estimated to be (7.1±0.6)%, while for 
Z° — > bb it is estimated to be (6.9±0.6)%. The combined signal expectation in OPAL is a factor of 1.1±0.1 times 
that of DELPHI, although the mass resolution in OPAL is 15 MeV and the background level is about 50% higher, 
when the b and c analyses are combined. 

The combined result of the bottom and charm analyses is shown in Figure ||, where events that pass either of the 
two analyses enter the histogram once. The Monte Carlo expectation for a signal the size of DELPHI'S observation 
is also shown. Limits have been obtained on the D*' + production rate by defining a mass window between 2.59 and 
2.67 GeV (±2cr), and comparing the observed candidate count against a background from a smooth parameterization 
shown also in Figure^. Systematic uncertainties, summarized in Table are incorporated in the limits. The 95% CL 
upper limit on the number of D*' + candidates in the combined analysis in the mass window is 32.8, with 443 events 
observed and 475 expected from the background fit. This can be converted into a limit on the production fraction: 

/(Z° -► D*' ± (2629)) x Br(D*'+ -> D* + 7r + vr-) < 2.1 x 10~ 3 (95%CL). (4) 
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FIG. 3. Reconstructed mass of D* /=t candidates in OPAL data (a), and in Monte Carlo simulation (b). The D l± mass has 
been constrained to the world average when reconstructing the D* ;± candidates (see text). The normalization of the signal in 
(b) is with respect to the DELPHI observation, and the width is taken to be 0.0 GeV. The solid error bars in (b) represent 
the true Monte Carlo statistical error, while the thin bars correspond to the errors of a sample the size of the data sample. A 
smooth parameterization of the background is shown for both the data and Monte Carlo, excluding the expected signal region 
indicated with arrows. 



No signal is seen separately in the tagged b or c samples, and the upper limits are 

/(c -> D*'+(2629)) x Br(D*'+ -> D* + 7r + 70 < 1.2 x 1CT 2 (95%CL) 

and 

/(b -► D*'+(2629)) x Br(D*'+ -► D* + vr+7r-) < 1.0 x 10~ 2 (95%CL). 



(5) 



(6) 



The analysis has been checked by reconstructing the orbitally excited D° and D2 , following a very similar selection 
procedure, except only one pion is added to the D*+ candidate. The rates and widths of these are consistent with 
Monte Carlo expectations. Combining the limits above with the measurements of these rates, which share common 
systematic uncertainties, OPAL sets a limit of R < 0.21 at the 95% CL, where R is denned in Equation |[ 



III. MEASUREMENT OF /(b -> A b ) x Br(A b -» AX) 

The measurement of f(b — > Ab) x Br(Ab — > AX) has been carried out mainly because of its interest in interpreting 
measurements of f(b — ► Ab) x Br(Ab — * A£~VX) to extract the quantity Rai = Br(Ab — > A£~77X)/Br(Ab — * AX), 
which is an approximate measure of semilcptonic branching fraction of the Ab- The average lifetime of the Ab has 
been measured to be approximately 80% that of the B mesons, while theoretical predictions can accommodate ratios 
between 90% and 100% jL5|. It is therefore important to check the shortened lifetime with a measurement of the 
scmileptonic branching ratio, which is expected to be proportional to the lifetime within the B hadron family provided 
that the leptonic decay currents are the same for the b baryons and the B mesons. A second reason for measuring 
this product branching ratio is that it allows one to extract Br(Ab — » AX). 

This analysis is fully described in Reference |l6[ ] and a summary is provided here. The procedure is to identify 
A ^p7r~ decays in b-tagged events, and then to use the momentum, p, and transverse momentum, p t of each A 
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TABLE II. Summary of the systematic uncertainties on the selection efficiencies for D*' + in the c and b samples. 



error source 


c sample 


relative contribution 
b sample 


combined sample 


relative errors on Monte Carlo efficiency 








Monte Carlo statistics 


9.3% 


9.1% 


6.2% 


detector resolution 


8.9% 


0.7% 


5.1% 


kaon dEdx 


U.o/o 


n no/ 
U.y/o 


U. / /o 


pion dEdx 


U.o/o 


U.o/o 


n ao/ 
U.o/o 


i j widen 


1 9<K 


ID. / /o 


o.D/o 


tt"q rrmont a i"ii~\n vni~\r\ oil n rr 
11 clJ^iIUJlllicLLlUIl IllUUdillg 


9 7% 


^ .o /o 


1 q% 


production rates in b/c 


none 


none 


0.3% 


total efficiency error 


13.2% 


19.3% 


9.1% 


relative errors on external branching ratios 








branching ratio D* + ->D°7T + , D° -^K~ty + 


3.1% 


3.1% 


3.1% 


error on r q q/r had 


4.5% 


0.6% 


none 


total branching ratio error 


5.5% 


3.1% 


3.1% 



with respect to its jet, to separate the contributions from b baryon decay, B meson decay, and A's produced by the 
fragmentation process. 

Events are tagged as Z° — > bb if at least one jet tags as a b using displaced secondary vertices as described in [ p"7| . 
The estimated overall efficiency is about 21% per jet and the purity is about 96% per event. The A selection requires 
that the reconstructed mass of the p7r~ combination be within 8 McV of toa, that it be farther than 6 MeV from 
the K° mass when interpreted as a 7r + 7r~ pair, that the flight distance be at least 8.0 cm from the interaction point 
when projected into the xy plane, that the momentum points back at the interaction point, that it have a momentum 
exceeding 5 GeV, and that it be within 0.2 radians of its jet axis. In addition, OPAL's dE /dx capability is used to 
require that the proton's candidate track's ionization is consistent with that expected from a proton. 

The Monte Carlo is used to provide distributions of p and pt for A's from each of the three sources, b baryon 
decay, B meson decay, and fragmentation. Other backgrounds contribute 3% of the sample, dominated by D + mesons 
tagging as B mesons because of their long lifetime. Because the D + does not have enough mass to decay into a A and 
another baryon, the main contribution comes from fragmentation A's in events with tagged D + mesons. These are 
incorporated into the fragmentation lambda portion of the fit. These distributions of p and pt are jointly fit for the 
fractions of A's in each of the three categories, and the results are shown in Figure [| and summarized in Table [II . The 
momentum is a powerful discriminator between b baryon decay and fragmentation, while the transverse momentum 
is effective for separating b baryon decays from B meson decays, although this separation is more difficult, which is 



reflected in the correlations shown in Table III . A summary of the systematic errors on the fit fractions is given in 
Table IV . Checks of systematic biases in the p and p t spectra were made by comparing events with A-lepton pairs 
and comparing the averages of p and pt in the right- and wrong-sign samples to obtain estimations for a pure sample 
of semileptonic A decays. 

The other important ingredient in the analysis is to understand the efficiency of the b-tag for jets containing a At,. 
This can be accomplished by counting events with one and two b-tagged jets while assuming the value of i?b and the 
efficiency to tag non-bb events and solving for eb, the b-tagging efficiency, using 



Rh^h + (1 — Rb)^udsc — fl 



and 



(7) 



Rh^h 



(1 - Rb)el dsc - h 



3 The p t of the A is measured with respect to the nearest jet axis, the calculation of which includes the A momentum. 



G 




p(A) GeV/c p t (A) GeV/c 

FIG. 4. Distributions of A p and p t in data (filled circles with error bars) and Monte Carlo (histogram). The three sources, 
fragmentation, B meson decay, and A b decay are normalized to their fitted fractions. 



TABLE III. Results of the fit for the fractional contributions to the A sample in b-tagged events. The correlations between 
the fit fractions are also shown. 



Source 


Fit fraction (%) 


A b -> AX 


Correlations 
frag^ AX 


B^ AX 


A b AX 


37.4 ± 5.3 


1.0 


-0.15 


-0.64 


frag^ AX 


37.1 ± 5.1 


-0.15 


1.0 


-0.61 


B^ AX 


25.5 ± 6.6 


-0.64 


-0.61 


1.0 



TABLE IV. Systematic uncertainties on the Ab fit fraction, fh b ■ 



Sources of Systematic Errors for /A b 


negative 


positive 




errors 


errors 


p(A) and p t (A) from B-^ AX 


-3.5% 


3.5% 


p(A) of Fragmentation A's 


-1.3% 


0.8% 


Pt (A) of Fragmentation 


-12.5% 


12.7% 


p(A) from A b -> AX 


-2.9% 


4.8% 


pt (A) from A b -> AX 


-16.6% 


19.8% 


Tracking Uncertainty 


-2.7% 


2.7% 


Total 


-21.5% 


24.4% 
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TABLE V. Quantities needed for evaluating the product branching ratio, and their impact on the uncertainty of the mea- 
surement. 

Uncertainty on 

/(b -> A b ) x Br(A b -> AX) 

N A 1582 ± 40(stat.) ±0.07(stat.) 

/ Ab 0.374 ±0.053(stat.)±°;°^,(syst.) 0.38(stat.)lo j?(syst.) 

R h 0.2169 ±0.0012(syst.) ±0.01(syst.) 

N mh 2323302 

e A 0.117 ±0.006(syst.) ±0.14(syst.) 

e b sig 0.294 ±0.015(syst.) ±0.14(syst.) 

Br(A -> p?r~ 0.639 ± 0.005(syst.) ±0.02(syst.) 



Large biases enter, however, from the facts that 1) the Ab has a shorter lifetime than the B mesons and therefore 
tags less often, 2) the A decay products are not included in the secondary vertex, leaving a lower average multiplicity 
and momentum of charged tracks which do participate in the secondary vertex, and 3) that if a high-momentum 
fragmentation A is present, the b hadron will have less momentum than average. These biases are accounted for by 
counting the number of selected A's in hemispheres with b-tags and hemispheres opposite b-tags. The tag efficiency 
corrections for each class of A are modeled in Monte Carlo, and the overall efficiency correction is scaled to match 
that of the data, obtained from the A-tag statistics. The scaled correction for the Ab — > AX component is then used 
as the correction, which amounts a factor of 0.49 ± 0.06, for an estimated event b-tag efficiency of (29.4 ± 1.5)%. 
The product branching ratio is then formed as follows: 



NAfA b 

2R h N mh e A e b slg Br(A -> P ti-)' 



/(b - A b ) x Br (A b - AX) = on Ar b AJ £, k — , (9) 



where is the number of A candidates in the final sample, N m h is the number of hadronic Z° decays in the final 
sample, e b sig is the b-tagging efficiency for events with Ab's in them. Table [v] gives a summary of the uncertainties 
on the measurement of the product branching ratio. Combining these factors with the appropriate systematic errors 
added in quadrature yields 

/(b -> A b ) x Br(A b -» AX) = (2.67 ± 0.38(stat.)±°;^( S yst.))%. (10) 

Combining this measurement with a previous measurement of the product branching ratio using a companion baryon 
technique jl^] yields a more precise estimation, as the overlap in the selected event sample is less than 20%. The 
combined measurement is 

/(b -> A b ) x Br(A b -► AX) = (3.50 ± 0.32(stat.) ± 0.35(syst.))%, (11) 

in agreement with a measurement from DELPHI jlj| of (2.2lJ |)%. Using the measured value of /(b — > Ab) in 
hadronic Z° decays of (10.1±|;?)% @, one obtains 

Br(A b -> AX) = (35±^)%. (12) 



IV. SUMMARY 



OPAL has produced a wide spectrum of heavy-flavor physics results in the past year. In addition to the measure- 
ments described here of the b baryon longitudinal polarization, the search for the radial excitation D* /+ , and the 
measurement of the product branching ratio /(b — * Ab) x Br(Ab — > AX), OPAL has completed its final i?b measure- 
ment |2lj, studied CP violation in the J/^-K^ channel p2|] , measured the B ± and B° lifetimes using a topological 
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tag |23], and measured the semileptonic branching fraction of inclusive b hadrons 24 1. Each analysis uses OPAL's 



full LEP1 data sample and are among the most precise published measurements of their respective quantities. 
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